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Abstract 



In terms of the s-wave phase shift of the two-body scattering at 
thermal wavelength, a systematic perturbative expansion of the Virial 
coefficients is developed for two-dimensional dilute system of bosons in 
its gaseous phase at low temperature. The thermodynamic functions 
are calculated to the second order of the expansion parameter. The 
observability of the universal low energy limit of the two dimensional 
phase shift with a quasi-two dimensional atomic gas in an anisotropic 
trap is discussed. 
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1. Introduction. 



Virial expansion in two dimensions 



The recent advent of laser cooling and atomic trapping techniques makes the physics 
of a dilute quantum gas experimentally accessible, which led to the observation of the Bose- 
Einstein condensation of the trapped metallic atoms in three dimensions [1]. Theoretical 
interests on these many-body quantum mechanical system has also revived since then. 
While the Virial expansion was investigated for a three dimensional dilute gas in both 
gaseous phase and condensate phase long ago [2], [3] and [4], a parallel formulation in 
two dimensions remains to be developed. There are many elegant works on 2D bosons 
concerning the quasi-Bose condensate near the absolute zero [5] [6] [7] [8] and [9]. The 
Virial expansion developed in this paper is complementary. As the quasi-two dimensional 
gas of trapped atoms is also experimentally feasible, the result reported in this paper may 
be brought to a direct comparison with the measurements 

As is well-known, a perturbative treatment of a dilute Bose gas in two dimensions 
suffer from two difficulties: 1) The scattering amplitude vanishes in the zero energy limit 
and the Born expansion breaks down for a large number of potentials [10] [11]. 2) The long 
range order parameter corresponding to the Bose condensate ceases to exist at nonzero 
temperatures because of the fluctuation of the condensate phase [12]. Both diflBculties 
stems from the two dimensional character of the density of states at low energies. We shall 
focus on the first difficulty in this paper. 

The Hamiltonian of a dilute system of interacting bosons is given by 

H ^Yli^P^ ~ + ^ 1] ,pi +Pi < 1 ^ |PiP2 > hl,^ hp^ hp^ ( 1 ) 

P Pl ,P2 ,P[ ,P2 

with bp, bp. the annihilation and creation operators of the bosons in their momentum space 
and jj, the chemical potential. We have chosen the unit of mass such that the mass of a 
boson is |. In the rest of this paper, two-body potential V is assumed to be isotropic, 
repulsive and of short range in coordinate space. 

The physics of the system depends on the relations among three length scales, the 
range of the interaction, tq, the average inter-particle distance and the thermal 

wavelength, 

A = (2) 

with (3 = {kBT)~^. The diluteness is measured by nrg << 1 and the classical limit cor- 
responds to 0. The Virial expansion we shall derive applies to the temperature 
such that ~ 1 when quantum coherence becomes significant. The effective expan- 
sion parameter is the s-wave scattering phase shift at the thermal wavelength, which is 
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^( in(A^r^)-i )• Therefore, the interaction corrections are far more significant than that of 
a three dimensional Bose gas at the same diluteness. 



In the next section, we shall introduce a renormalized two body potential 
pansion parameter in the dilute limit, which set up the systematics of the perturbation. 
The corresponding Virial expansion to the second order of the renormalized potential will 
be developed in the section III with typical thermodynamic quantities calculated in the 
section IV to the same order. The convergence of the Virial expansion and its relation to 
the quasi Bose-Einstein condensation will be discussed in the final section. 



2. The Renormalized Potential. 



The requirement of renormalizing the interaction potential is not unfamiliar in three 
dimensions with a hard sphere potential or more realistic Lennard-Jones potential, for 
which a straight forward perturbation series breaks down because of the singular behavior 
of the potential at short distance. The natural choice for the renormalized potential is the 
exact two-body scattering amplitude at zero energy (scattering length). What is lacking 
in two dimensions is such a natural choice, since the scattering amplitude vanishes at 
zero energy. This is analogous to perturbation theory of the quantum chromodynamics, 
where the infrared slavery and the asymptotic freedom deprive us of a natural scale of 
the ultraviolet renormalization. A running coupling constant defined at relevant energy 
scale has to be introduced as the expansion parameter. For a two dimensional Bose gas, 
we need also to introduce a running coupling constant, which turns out to be the s-wave 
phase shift at the thermal wavelength of the bosons. 



To replace < p*^P2l^|PiP2 > with appropriate renormalized quantity at low energies, 
we focus on the two body sector of the Hamiltonian (1) and define 



+ e 



(3) 



The hermitian operator V can be formally expanded according to the power of the potential 
V and vice versa. To the leading order V = V and it is straight forward to show that 



V = V + 



1 



(4) 



Sandwiching V between two-body states and completing the integral over r, we find that 



<PlP2\V\PlP2 >=< PiP2\^\PlP2 > + 



pi+pi-kf-ki pf+pi'-kf-ki 
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X < ^^^^\V\kik2 >< kik2\V\pip2 > +0{V^). (5) 

As we shall see, the matrix element < p'iP'2\'^\PiP'i > vanishes in the limit /? — > oo and is 
a proper choice of the renormalized potential at low temperature. Furthermore, replacing 
V by its formal expansion in V, the systematics to all orders is restored. In terms of the 
two-body matrix element of V, the Hamiltonian of the system, (1) becomes 

H = J2^P^ ~ l^^^ + ^ Yl ^P^+P^ 'P'l +P2 < -^1^2 \V\piP2> bi,^ 6t, bp^ hp^ 



2n 



+ ^ +P2 ,p[ +^^2 < 1 V - 1/ \p1p2 > 6t, bp^ bp^, (6) 

where the last term is understood as an power scries of V starting from the order V^, 
like the renormalization counter terms in relativistic field theories. A similar method of 
renormalization was developed in the context of a lattice gas with on-site exclusion. [13] 

The matrix element (3) is related to the binary kernel of Lee and Yang [3] through 

< P1P2W1P2 >=-l[< Al^2\BmPiP2 > e^(P?+^^)+ < PiP2\B{PMp'^ >* ef'^P'"+p"\ 

(7) 

and the Virial expansion developed here for a 2D Bose gas is parallel to that of Lee and 
Yang for a 3D hard sphere gas. 



In terms of the total momentum P = pi +P2, P' = p[ and the relative momentum 
P= -P2), ^ = Kp^i -^2)1 we have 

< P1P2MP1P2 >= Ih'M^^'^"' < 0\Ve-^''\p> +e^^p" < 0\e-^''V\p>], (8) 

where 

h = -2Vl + V{r) (9) 

is the Hamiltonian of a potential scattering problem with r the relative coordinate. In- 
serting the complete set of eigenstates of h, i.e., h\n >= En\n >, we have 

< p'lFe-^'^lp >= ^e-^^" < ^\V\n >< n\p> (10) 

n 

and < '^\e~^^V\p >=< p\Ve~^'^\'^ >*. For an isotropic potential, \n > is specified by 
the azimuthal quantum number, m, and the radial momentum, k, with E = 2k^. The 
corresponding wave function 

< f\k,m>= Nkm^e''^'^um{k\r) (11) 
V 27r 
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where the radial wave function Um{k\r) approaches asymptotically 



Umir\k) c:^ ^ ^ cos kr - ^ - J + 5mik)^ (12) 

for A;r >> 1 with Sm{k) the phase shift and A^^^ is a normalization constant. The partial- 
wave expansion of (10) reads 

1 2 °° 

< ^ I Ve-^^ \p >= -Vo{p',p) + -J2 ip', P) cos m$, (13) 

m=l 

where 

/■OO nOO 

K.(p',p) = 27rVAr2^e-^^^' / dr'r' Jm{p'r')V{r')um{k\r') / drrUm{k\r)Jm{pr) 
k -^0 

(14) 

with $ the angle between p and p', and Jm{z) the Bessel function. The infinite volume 
limit of the integral over r and the sum over k in (14) have to be evaluated carefully. We 
defer the details to the appendix A and quote only the result here, 

Vm{p',p) = -4e-2^^Y^)"^sin25^(p) + 7rP H dkke-'^^' ''^^^':^^''^^^'^^ (15) 

\p / Jq k —p 

with 

Vmip,k)^ drrJmipr)V{r)um{k\r), (16) 
Jo 

4 

Vm ip, P) = sin 5m ip) ■ (17) 

TT 

and V the principal value of the integral. For a short-range potential, the wave function 
UTn{k\r) normalized according to (12) takes the form 

^m{k\r) c '^^^fmir) (18) 

in the limit kr 0, with fm{f) independent of k. The phase shifts display the following 
low energy behavior 

and 

Sm{k) - A:^-, (20) 

where a is the s-wave scattering length in two dimensions. For a hard sphere potential, 
a = \e^rQ with tq the radius of the sphere and 7 the Euler constant. But the scattering 
length and the range of the potential may not be comparable in general. The relevant 
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length for the Virial expansion developed in this paper is the scattering length a. It 
follows from (17) and (18) that 

Vm{p,k) = -^{^^sin5m{k) (21) 

for k and p both small, and the low energy and low temperature approximation of (15) 
reads 

Vm{p\p) ^ -4e-'^^U^-T s\n2Smip) + -p'^p^^V H dkk-'-^+'e''^'^''^^^ 

\p / TT Jq k^ — p^ 

(22) 

The eq. (19) was proved rigorously by Chan, Khuri, Martin and Wu [11] for a general 
class of potentials that fall off faster than — for r ^ oo . They also proved that for the 
same class of potentials, the correction to the corresponding function /o(r) of is of the order 
of /c^. As this involves only the long wavelength limit of the scattering, their conclusion 
can also generalized to the repulsive potential that becomes singular as r — > 0. Therefore, 
all the logarithmic dependence on the range of the interaction of. In a, is absorbed in the s- 
wave phase shift through (22). This way the s-wave channel dominant over all other partial 
wave channels to all orders of a low energy expansion, different from three dimensions. The 
universality found in [11] is highlighted in the low temperature thermodynamics of the 2D 
system. In what follows, we shall suppress the subscript "0" for the s-wave and introduce 
a running coupling constant at the scale of the thermal wavelength, 

«=— 1^^ , (23) 



we have a ~ ;^(5( ^^^^ ) and the s-wave phase shift becomes 

5 = -TTCK + TTCt^ ( In h7)+... (24) 

V 27r / 

AS we shall see, the choice of the constant pertaining to the logarithm of (23) is to make 
the second Virial coefficient free from 0{a'^) corrections. It follows from (10), (15) and 
(24) that the renormalized potential at low energies reads 



dp p, r / rm X \ Q^ii{p-ii-) 

< A^2MpiP2 >= [Sttck - Sttck^ (^In + 7 j + IQTra^V J dkk ^2 _ ^2 + 

and this expression will be applied extensively in the subsequent sections. 



The dependence of the dimensionless running coupling constant (23) on T and a is 
completely analogous to that of a relativistic field theory of zero masses in 3D with T 
corresponding to the renormalization energy scale and a the ultraviolet cutoff. 
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3. Virial Expansion. 



Virial expansion in two dimensions 



The thermodynamics of a uniform gas is determined completely by its equation of 
state, usually expressed in the form of the Virial expansion, 

oo 

^ 1=1 

and 

with n the number density and z — e^^ the fugacity. The l-th Virial coefficient, bi is 
determined by the quantum mechanics of I particles, the exact solution of which is in 
general unavailable for I > 2. For an ideal Bose gas in two dimensions, bi = On 
writing 

^1 = ^ + ^1^ (28) 

we have 

^ \92{z)+T (29) 



with g2{z) = Yl'ili F ^ = Yl'ili M-^'- The perturbative expansion of F is represented 
by thermal diagrams. 



The thermal diagrams of F to the third order in a is shown in Fig. 1, where a solid 
line represents a boson propagator, -p^+n with tOn the Matsubara energy of the boson, 
a solid circle vertex is associated to the factor < p'iP'2\V\pip2 > with pi, p2 (p\, p'2) the 
incoming (outgoing) momenta and an open circle vertex denotes < ^1^2!^ ~ V\pip2 >, 
the analog of the rcnormalization counter term in a relativistic field theory. In what 
follows, we shall calculate the F to the order a^, which is one order beyond the mean field 
approximation. 



On writing 

F = F« + F6 + Fe + Fd + ... (30) 

with Fa, Fb, Fc and F^ standing for the contribution from the first four of the diagrams in 
Fig. 1 in a sequential order, we have 

= --^ > Y- 2^ — \r- \ — V / / 7T~\2 <PlP2 VpiP2 > n{pi)n{p2) 

/32 ^ [i^^-p^^ n){i^^-p^^ J (27r)2 7 (27r)2 



72- 



/Oj ~* f u ~* 
(2^ J (2^ <PiP2MpiP2 > n{pi)n{p2) (31) 
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with 



n{p) = 



ze 



1 - ze-^P'' 



z e 



-Wp^ 



(32) 



and z = e^^ the fugacity. In accordance with the expansion (25), we have 

r — r' -I- r" 

with 



r'^ = -8na(3 



d pi f d p2 



■n{pi)n{p2) 



1 /I P'^\ 



and 



(27r)2 7 (27r)2 

—T7^n{pi)n{p2)V / dkk^ — 



d'^Pi r d'^p2 

(2^7 (2^'"^"^^'"^"^^' 7o ^"""p^-k^ 
Using Taylor expansion of (32) and integrating over pi and p2, we end up with 



(33) 
(34) 

(35) 



A2 



In^ 1 . ^ . -.^2, 



1-^ 



+ Q;i:>(^) + 0(q;^) 



(36) 



with 



^^+* , 2rs 



In 



r,s=l 



rs r + s 



Similarly, we have 

1287r2a2 f d'^pi f d^p2 f d'^pa 



f d'pi f d'p2 f d'ps ^ 
J (27r)2 7 (27r)2 7 (27r)2 2^ 



{iuji -pl+ iJ,)^{ioJ2 -pI + IJ'){i(^3 - + A*) 



8a^ z , 2 1 



In^ 



Finally 



A2 (1-^)2 1-^ 
(37) 



2 2o f d^Pl f d P2 f f d^p^ 2e2 



X 



(27r)2 7 {2t^Y J {2TTy J (27r)2 



1 



■ (27r) 5 (pi + P2 - P3 - pk) 



0)1,0)21^3 



(iwi - + /u) (iw2 - + /") (^t^3 - pI + A*) W'^i + Ci;2 - ^3) - P? + A*] 



327r2Q;2/3 f d'^pi f d^p2 f d^ps f d^p^ 



72- 



72-. 



(27r)2 ,/ (27r)2 ,/ (27r)2 j (27r)2 



(277)2^2 ^p-^ -P3 -P4) 



P?+P2 -P3 



n{pi)n{p2)n{p3)n{p4), 



(38) 
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which represents a genuine three-body scattering, and 



= -64.V^y ^P-)^^>-' -p, -p.) 

e/3(p?+P2-P3-pI) _ 1 

^ „2 , „2_„2_„2 MPiMP2), (39) 
Pi + P2 P4 

Things are greatly simphfied to the order Oip?^ by forming the following combination 

r', = r;' + re + rrf 



(40) 

This expression is then simplified in four steps: 

1) Expanding n(pi) according to the power of z\ 

2) Using the property of the principal value, 

1 /"^^ 
= Re = = Im / f]rp^(pi+P2-pl-pl+iO'^)^ (41) 



1 



pI + P2 -pI-pI 



and the Fourier transformation of the delta function 

72;? 



5\pi +P2-Ps-Pd- J ^e^(P^+P^-P^-P^yP; (42) 



3) Carrying out the Gauss integration over jf s and then the Gauss integration over p 
for each term of the power series in z; 

4) Carrying out the elementary integral over x. The final result reads 

4^2 

r'e = -^Hz) (43) 

with 

r,t^l + 0(« + \/(^ + t){s + t)- VFi' 
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Collecting above results, we obtain the Virial expansion of a dilute Bose gas to the 
second order of the interaction, 



P 1 / N 



2a, 2 



i^ln^ 



1 



with 



1 



0(z) ^ F{z) + -D{z), 



and the corresponding number density is given by 

Aa z 



n 



1 1 1 



A2 1 



1 d 
In- h 



1 



z 



A^ din 2 



i^ln^ 



+ 20(z) +0{a^), (44) 



(45) 



+ 20(^) +0(a3). (46) 



We notice that the order ct^ term start with the third power of z, which may be viewed as a 
criteria to fix the constant part pertaining the logarithm of the running coupling constant 
(23). Furthermore, the coefficient of agrees with the result obtained with the classical 
formula of Beth and Uhlenbeck [14]. The asymptotic behavior of D{z) and F{z) as 2; ^ 1~ 
is analyzed in the appendix B. 



4. Thermodynamical Functions. 



Among experimental observables of a two dimensional Bose gas are the homogeneous 
thermodynamical functions, which follow readily from the Virial expansion of the last 
section. For the sake of clarity of notations, all thermodynamic functions in their ideal gas 
limit will carry the superscript '(0)'. 

Inverting (46), an expression of the fugacity z in terms of the density n is obtained to 
the order a^, 

Aa^ d(f) 



\nz = Inz^^^ +4q;^ - 



— 1 din 2; 



s(o) 



(47) 



with ^ = nX^ and z^^^ = 1 — e ^, the fugacity of an ideal Bose gas in two dimensions. On 
substituting (47) back to (44), we obtain the equation of state to the order a^. 



p = + 8W + U(l - e-^) - 

A^ L — 1 a m 2 



2=l-e-« 



(48) 



where 



P 



(0) 



47r 



47r 



(-) 



m— 1 



(2m + 1) 



^^^m+ij (49) 



is the equation of state for an ideal gas with the m-th BernouUian number. 
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The Helmholtz free energy per unit volume is obtained through the formula 



/ = fj^n — p = f^^^ + Snan^ — 



A4 



-0(1 -e-«) 



The entropy per unit volume is given by 

167ra^ ^ 

= s^"^ - 



e - 20(1 - e-^) + 



- 1 din; 



z=l-e-S. 



(50) 



(51) 



and the correction to the ideal limit is of the a^. It follows from (50) and (51) that the 
internal energy per unit volume is 



— 1 d In 2; 



2=l-e-f 



(52) 



The specific heat at a constant volume (area) reads 



cv 



,(0) 



A2 



^(2e^-2 + eeS) d<t> 



(e^-l)2 dlnz {e^ -ly d{\YvzY 



2=i-e-e 
(53) 

and the leading order contribution of the interaction is proportional to . The isothermal 
compressibility, kt can be calculated readily from (48) with the results 



= _^(^\ ^ ^(0) + 16W + 



d(j) 



'KdQJT " ' (e^ - 1)2 r dln^ d{hizY 

Employing the thermodynamic relationship 



j=l-e-« 



Kg - Kt = 



'T/dp 



) \-(-) 



Cv 



we find the adiabatic compressibility 



dp \ 



K5 = — nf-^) = kI^'* + 167rQ!n^ + 167rQ:^n^ + 
\dQ.J s 



2„2 



327ra2 



4>{z)-- 



— 1 d\n.z 



z=l-e-i 



which can be directly measured through the sound speed in the Bose gas, 



V = 



2ns 



n 



(54) 



(55) 



, (56) 



(57) 



5. Concluding Remarks. 



In previous sections, we have developed a diagrammatic approach to the Virial ex- 
pansion of a dilute Bose gas with a repulsive interaction of a short range, valid in the 
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region where nA^ ~ 1 and A >> a. Before concluding the paper, we shall remark on its 
implication on the quasi Bose-Einstein condensation and its applicability to the realistic 
system of an quasi-2D gas of alkaline atoms in a strongly anisotropic trap. 

All coefficients of the Virial expansion (44) are singular at 2; = 1, and the singularity 
gets stronger at higher orders. A simple power counting argument shows that the coefficient 
of diverges as (1 — zY~^ as 2; — > 1~, up to some powers of In y^^, as is bore out in 
the explicit form of (44) to (See appendix B for the singularities of D{z) and F{z)). 
Therefore the reliability of the expansion requires that 

a«l-z (58) 

Using the formula (23) for a and the ideal gas limit of 2;, the condition (58) becomes 

T»y^. (59) 

which is consistent with the formula of the transition temperature to the superfluid phase, 
obtained in [6] and [8], under a mean field approximation. It will be interesting to see 
whether the resummation of diagrams developed in [15] for a 3D hard sphere Bose gas can 
be applied here to extract the ground state energy beyond the mean-field approximation 
and to compare with the results in [7] and [9]. At this stage, we only remark on that 
the order a correction to the internal energy is twice of the rigorous result of Lieb and 
Yngvason [16] for the ground state energy for the thermal wavelength comparable with the 
inter-boson distance. The factor two comes from the exchange energy which is absent at 
zero temperature because of the Bose condensate. For a point-like repulsion corresponding 
to the leading order term of (25), the exchange energy is equal to the classical energy 
and thereby doubles the interaction energy under the mean-field approximation. A similar 
effect was observed in three dimensions [15] t 

The Virial expansion in two dimensions is characterized by the logarithm in the denom- 
inator of the running coupling constant (23) with a universal coefficient. The observation 
of such a universality demands the logarithm to be large enough a large enough such that 



t Another way to understand the factor two is to consider the expectation value of the 
potential energy of a point-like repulsion with respect to a state of A/'(= nO) free bosons, i. 
^- E =< \^T.p,,p2,p[,p'^ ^Pi+P2,p',+P'2^f^^f^hih2 \ > with I > a product of single particle 
states. The summation can be broken into three nonvanishing terms: E = -^ly^fj ^ < 

KAM^P2\ > < I^J-,^kViV2l > -Ep < \b^^^^pbp\ >]• In the limit of an 

infinite volume, — > 00, each of the first and second terms gives rise to jV^. If all bosons 
occupy one momentum level (Bose condensate), the third term contributes — jV^, making 
the sum ^ = ^Vn?. If none of the levels are macroscopically occupied, the third term is 
0{J\f) and will not contribute to the infinite volume limit. We have then ^ = Vn^. 
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the first few terms of the expansion represent a reasonable approximation. Experimentally, 
a quasi-two dimensional gas can be implemented with a strongly anisotropic trap, which 
can be modeled as a three dimensional gas in a narrow harmonic well in one direction, 
referred to as the trapped direction. An analytical solution to the two-body scattering for 
large extension of the single particle wave function in the trapped direction, Z, in compari- 
son with the range of the inter-particle interaction (which remains three-dimensional) was 
obtained in [17], and an effective phase shift of the two dimensional s-wave component can 
be extracted when the de Broglie wavelength in the trapped plane becomes much longer 
than /, 

5efr. = , (60) 

2 In — 

where k is the relative 2D momentum and tteff. = y^^^xp ^ — a/^^^ with a the 3D 

scattering length. The authors of [17] also found numerically that the approximation (60) 
works well even when I ~ o. The typical inter-particle distance in an atomic trap is about 
lO'^A and the typical 3D scattering length for alkaline atoms is of the order of lOOA. For 
the thermal wavelength equal to the inter-particle distance and the trapped dimension 
equal to a (which is technically feasible), we have tteff. — 5lA and then a ~ 0.12, according 
to (23) with a there replaced by Oeff.. The universal 2D logarithm could be quite significant 
to observations. 

Acknowledgment I would like to thank Professor N. Khuri for valuable discussions. 
This work is supported in part by the US Department of Energy under grants DE-FG02- 
91ER40651-TASKB. 

Appendix A 

To regularize the infinite volume limit of Vm{p' ,p) defined in (13), we follow the 
usual practice by restricting the relative coordinate f within a large circle of radius R 
and imposing the Dirichlet boundary conditions for the wave functions with and without 
potential V, i.e. 

Ura{k\R) = (Al) 

and 

ut\p\R) = Q, (A2) 
where Um{k\r) is the solution of the radial Schroedinger equation, regular at r = 0, 

- rTri'd^) + " V{r)]uUk\r) = 2/c\^(/c|r), (A3) 

and Um\p\r) = Jm{w) satisfies the radial Schroedinger equation of a free particle, i.e. 

2 d / d\ 2m 



r dr 
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For pr » 1, the asymptotic behavior of tt^^(p|r) is 

^m\p\r) - cos (pr - ^ - 9 (^^) 

and that of Um{k\r) for kr » 1 is given by (12). It follows from (Al), (A2), (A5) and 
(12) that 

P = Pn = y' + j + i)^ (^6) 

and 

k^kn^Pn ^ (-47) 

for large R. 



Multiplying (A3) by Jm{pr) and (A2) by Um{k\r), subtracting the result and using 
the boundary condition (Al) and (A2), we find 

r drrJm{pr)um{k\r) = J^^^^'^^, , (AS) 
Jo 2{k^-p^) 

where 

Vmip,k)= drrJm{pr)V{r)um{k\r). (A9) 

JO 

and is well behaved in the limit R^ oo. The eq. (14) becomes then 

k " ~p 

with A^fcm the normalization constant such that 

NL / drruUk\r) = 1. (All) 
Using the asymptotic behavior (12), we find Nk^n — A/ ^ ^r large R. 



To take the limit R ^ oo of the sum (AlO), we need to isolate out the /c's which 
are sufficiently close to p. For this purpose, we introduce a subset of /c's, B, such that all 
k e B satisfy the condition that \k — p\ < ^ for N » 1. We further specify the order 
of the limit such that i? — > oo first and then A?" — > oo. Consequently, the summation in 
(AlO) is divided into two parts, 

Vm{p',p) = V<{p',p) + V>{p',p). {A12) 
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with V^{p',p) including only the fc's within B and V^{p',p) all others. For V^{p',p), all 
the k^s except those in the denominator can all be set to p, and we obtain that 

TT^ 2 ^ 1 Tt"^ 2 

Vmip'^P) - Vm{p' ,p)Vm{P,P) —. = ^6"^^^ Vm{p' ,p)Vm{p,p) COt 5 m{p) ■ 

(A13) 

The summation Vm{p' iP)-, however can be replaced simply by the principal value of an 
integral, i.e. 

Following the same steps that leads to (A8), we derive that 

4 

Vm {p, p) = sin 5m {p) ■ {Alh) 

TT 



Appendix B 

In this appendix, we shall determine the singular behavior of the functions D[z) and 
F{z) in the Virial expansion (44) as 2; ^ 1~. 

On writing D{z) = Di{z) + D2{z) with 

00 ^r+s 



Di(z)= y ln2rs (SI) 



and 

we have 
where 

and 



, rs 

r,s=l 



°o ^r+s 



D2{z) = -y ln(r + s), (S2) 

f ^ ret 



rs 

r.s—l 



Di (z) = 2Si {z)S2 (z) + Sf (z) In 2, (S3) 



n=l 



00 ^ 11 

^2(^)-E^i^^^2^^' w ^^^^ 

n=l 
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as 2; — > 1~. To estimate D2{z), we rewrite it as 

00 

D2iz)=J2CNz'' iB6) 

N=l 

where 

A''—! 

Cat = lnAr V ^ = A IniV V - ~ |- ^ (B7) 

m+n=N n=l 

for iV » 1. It follows then that 

°^ N 9 1 

' ^ N 3 \-z ^ ' 

N=l 

Combining (Bl) and (B2) with the aid of (B3-5) and (B8), we obtain the asymptotic 
formula 

D{z)c^\ln^-^. (39) 
o 1 — z 

The asymptotic behavior of F(z) can be determined by its integral representation, 
read off from (40) and (43) 

=-128. /3 Pj p-jj(2.) 6 (p.+p.-P3-P4) ^.^^._^._^. 

= ^ J J2^^J J2^^J J2^j W?^^"^ iP^+P^-Ps-Pd 

n{pi)n{p2)[n{p3) + n{p^)] - n{ps)n{p^)[n{pi) +n{p2)] , . 

Pi+pi-Ps- pi 



On writing z — e the leading singularity of F{z) in the limit z ^ 1 can be extracted 
from 

then 



from the integration domain where pj < -A= for j = 1, 2, 3, 4 with e << 77 << 1. We have 



Fiz\ ~ 327r3 { TT (27r)^<^^(fci + fc2-fe3-fe4) 

da;a;K^(a;) (511) 



'3 

Owl 

e 

with ^"0(2;) the modified Bessel function of the second kind 
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